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We find numerical solutions of Einstein equations and scalar field equation for a global defect in 
higher dimensional spacetimes (> 6 ). We examine in detail the relation among the expansion rate 
H and the symmetry-breaking scale r/ and the number of extra dimensions n for these solutions. We 
find that even if the extra dimensions do not have a cigar geometry, the expansion rate H grows as 
r 7 increases, which is opposite to what is needed for the recently proposed mechanism for solving the 
cosmological constant problem. We also find that the expansion rate H decreases as n increases. 
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I. INTRODUCTION 


Recent observations suggest that majority of the Uni¬ 
verse is the unknown: 27% is matter (and only 4% is 
ordinary matter) and the rest is the cosmological con¬ 
stant (or dark energy) Hi Since the cosmological con¬ 
stant is interpreted as a vacuum energy these ob¬ 
servations indicate that the energy density of the vac¬ 
uum energy is p\ « (10“ 3 eV) 4 . On the other hand, a 
field theoretical estimate of zero point energy of the vac¬ 
uum yields the Planck energy « 10 18 GeV. Therefore, the 
mismatch between observations and the theory is huge: 
(10 _3 eV) 4 /M 4 ~ 10 -120 . This mismatch is called the 
cosmological constant problem (CCP) Ell and still re¬ 
mains to be solved. 

Dvali et aLQ has suggested the mechanism of diluting 
the cosmological constant by using brane world models 
with codimension greater than 2. In this mechanism the 
observed effective cosmological constant related to the 
expansion rate becomes small because whose energy is 
consumed to bend the bulk space even if the bare vacuum 
energy (brane tension) is as large as the Planck scale. 
The success of the mechanism rests on the conjectured 
relation in higher dimension 


H ~ M* 



l/(n-2) 


(1) 


where M* is the (4 + rc)-dimensional Planck mass and 
PA 4 is the four dimensional (bare) vacuum energy (brane 
tension). If the number of extra dimensions n is greater 
than 2, the expansion rate H is inversely proportional to 
the brane tension pA 4 - Then the smallness of the cosmo¬ 
logical constant could be explained by the largeness of 
the brane tension. 

Cho and Vilenkin (extending the analytic solutions in 
0) have recently constructed numerical solutions of a 
global defect in seven (n = 3) dimensional spacetime Q. 
They obtained numerical solutions of an inflating global 
defect if the symmetric breaking scale is greater than the 
higher dimensional Planck scale. Then the inflation rate 
is found to grow almost linearly as the brane tension is 


increased, which is opposite to what is needed to solve 
the cosmological constant problem 0 

The main purpose of this paper is to extend Cho and 
Vilenkin’s model to the arbitrary (4 + n)-dimensional 
one and find numerical solutions without “cigar ansatz” 
and to examine the relation among the number of extra¬ 
dimensions, the energy density and the expansion rate. 

We shall find that the conjectured relation <0 does not 
hold even by using our new bulk solutions in other extra 
dimensions n(> 2). However we find that the brane’s 
expansion rate is a monotonically decreasing function of 
the number of the extra dimensions n and can vanish at 
the specific dimension. 

In the following section, we introduce the model in Sec. 
Ini Results of the numerical integration are shown in Sec. 
mu Finally, we summarize our results in Sec. EH Some 
numerical details are given in Appendix El Cigar type 
solutions are reexamined in Appendix iBl 


II. MODEL 
A. Space-time structure 

In our model, the brane is assumed a 4-dimensional 
de-Sitter apace dS 4 , and the extra space is a spherically 
symmetric n-dimensional space R x S” _1 . Where the 
number of extra dimensions n is equal to or greater than 
2. The entire manifold is wrapped product of both spaces 
R x S" _1 x dS 4 , whose metric is 

ds 2 = dr 2 + C{r) 2 r 2 dUl 2 n _ 1 

+ B(r) 2 + e 2Ht ^ dx i2 ^j . (2) 

Here the coordinate of the brane is ( t , x 1 , x 2 , x 3 ) and H is 
the positive constant expansion rate. The extra space’s 
coordinate is (r, 0 1 ,..., 0 n - 1 ) and dfl „_\ is the metric of 
an n—1 dimensional sphere S” -1 . C\r)r and B{r ) are the 
radius of the extra space and the warp factor depending 
on r only. We adopt the Einstein-Hilbert action for the 
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space-time dynamics such that 

Se-h = -^ J d A+n x^gU, (3) 

where n 2 = 1 /M 2+n with M* being the (4 + n)- 
dimensional Planck mass. 
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B. Energy Momentum Tensor 


The global defect in the n-dimensional spherically sym¬ 
metric space is introduced to construct the brane, which 
is described by a multiplet of the scalar fields </> 1 with the 
action, 


S# = / d 4+ "a;v /3 5 


-hd A 4>'d A (l)i - v{4>) 


( 4 ) 


where capital letters (A,...) and small letters (*,...) run 
from 1 to 4 + n and from 1 to n respectively. Because we 
consider spherically symmetric solutions only, the scalar 
multiplet has been assumed to have a hedgehog config¬ 
uration, (j )* = </>(r)£ l / r - Here <j>(r) depends only on the 
radius coordinate r and £* represent for the Cartesian co¬ 
ordinates of the extra space. The potential of the scalar 
field V(<f>) has minimum at \(f> l \ = <f> = rj such that 
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Here G)r — 12 H 2 represents for the 4-dimensional Ricci 
scalar depending on the expansion rate of the brane. The 
prime denotes the differentiation with respect to r. The 
equation of motion of the scalar field is 


Vtf) = \{<t> 2 -V 2 ) 2 - ( 5 ) 

The energy density due to the scalar field may be re¬ 
garded as the brane tension. 


<P" + (n- 



B' 

~)~B 



- (™-l)^-A#/> 2 -77 2 ) = 0. 


( 9 ) 


Eq.© imposes the constraint when solving eq. @, @ 
and © as the second-order differential equations for B , C 
and 4>. 


C. Basic Equations 


III. NUMERICAL SOLUTIONS 


The Einstein equations and the equation of motion of 
the scalar field are derived from the action, Eq. © and 
Eq.©. The Einstein equations are 


nv — __ 

CT 11. - 


1 ^R 


4 B 2 
+ 3(n - 1) 

+ 0-4) 

+ n(n — 1) 


B" 

i 3 

+ 3 | 

' ° B 

f B ' + 

B'C 1 

l Br + 

BC 


C" (n —2)(n —1) (C' 

~C + 2 V’G 

C (n — 2)(n — 1) (1 

Cr 2 \r 2 ( 

1)</> 2 _A 2 2 2 

2 2 C 2 r 2 4 ' > 


( 6 ) 


We have solved the Eq. ©, ©, © and © nu¬ 
merically with the initial conditions 13(0) = C(0) = 1, 
B'( 0) = C"(0) = 0 and <j){ 0) = 0. These equations have a 
set of three parameters (n, kt], (k/A 1 / 2 )#). It is found by 
the numerical integration that the proper relation among 
them is obtained under the condition that the point of a 
singularity becomes as far as possible. We assume that 
only particular combination of these parameters gives 
a nondiverging solution. We shall call such parameters 
eigenvalues and regard regular solutions with eigenvalues 
as physical solutions. Solutions obtained from param¬ 
eters deviated from the eigenvalues have divergence in 
B or C. Similar situations are considered in [§| for the 
cigar ansatz which is detailed in Appendix ET We shall 
find yet another numerical solutions. The set of eigen¬ 
values forms the surface in the 3-dimensional parameter 
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space, whose shape also will be studied in the following. 
The discussion of technical details for finding numerical 
solutions is given in Appendix 0 


A. Asymptotic Solutions 

For H = 0, we can find an asymptotic solution analyt¬ 
ically which is obtained by solving Eq. ©, 0 , @ and 
0 at large r. If n > 3, the solution is 

<p{oo)=rj, ( 10 ) 


B 2 (o o) = constant, (11) 


C 2 (oo) = 1 - 


(^) 2 
n- 2 ’ 


( 12 ) 


where (nr]) 2 < n — 2. From Eq. (11121) . the sphere S n 1 
has a solid angle deficit such that 


27 t "/ 2 ( act ?) 2 

= r(n/2) ' n- 2’ 


(13) 


where T is the gamma function. As nr] approaches 
\Jn — 2 , the deficit angle consumes the entire area. If 
n — 2, C{ oo) can take an arbitrary constant. 


B. (kti) 2 < n — 2 Case 

For ( nr] ) 2 < n — 2, a non-singular solution exists when 
the brane is not expanding, H = 0. Then, the solution 
takes an asymptotic form given in the last subsection. 

We have solved the Einstein equations and the equa¬ 
tion of motion of <f> numerically in the range of [ 0 , r max ]. 
Here r max should be taken to be sufficiently large so that 
<j> takes a constant value given in Eq. m- The details of 
the method of numerical integration is given in Appendix 

El 

As an example, Fig. Q shows a solution with the pa¬ 
rameter (n, nr], (k/A 1 // 2 )iJ) = (3, 0.65, 0). The scalar field 
</> approaches 77 rapidly, which makes the defect core with 
the energy density {nr]) 4 approximately, and B(r) and 
C(r) approach toward constants as Eq. Ill Oil . (tTTI and 

C2». 


C. (kti) 2 > n — 2 Case 

For (kt]) 2 > n — 2, arbitrary H including H = 0 leads 
to a divergence in B or C at finite distance from the 
origin and the singularity is formed. We call this point 
r S in g . At the specific value of H } the distance to the 
singularity becomes as far as possible and the divergence 



FIG. 1: This graph shows a solution with the eigenvalue 
(n, kt], {k/X 1/2 )H) = (3, 0.65, 0). 



FIG. 2: This graph shows the solution with the eigenvalue 
(n, nr/, {k/X 1/2 )H) = (3,1.09, 0.003786056). B approaches to¬ 
ward 0 at finite rf. 


vanishes. We call this point rf. We note that a fine- 
tuning of the parameters is required to find this local 
peak of the distance. The details of the numerical method 
is given in Appendix fXl 

As an example, the numerical solution with the 
eigenvalue (n, kt], (k/A 1 // 2 ).H) = (3,1.09,0.003786056) is 
shown in Fig. [3 It is noticed that the B(r) vanishes 
at finite r but C(r) does not diverge. The solutions for 
other dimensions can be obtained and are shown in Fig. 

m 

Fig. El shows the relations between 77 and H with n 
fixed at some values. Each line approaches the point 
(kt],H) = (y/n — 2, 0). We find that H grows as 77 is 
increased. This tendency is similar to the Friedmann 
equation but opposite to the conjectured relation Eq. 0 . 

Fig. 0 shows the relations between n and H with ij 
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FIG. 3: (a), (b), (c) and (d) are for the eigenvalues 

{n,Kr 1 ,{K/X 1/2 )H) = (2,0.65,0.003467), (4,1.46,0.004330), 
(5,1.765,0.004535) and (6,2.025,0.004119) respectively. 



FIG. 4: The relation between p and H with n fixed respec¬ 
tively. In the region where H is small, we cannot find eigen¬ 
value, because the point of singularity is very far from the 
origin. H = 0 seems to be established at np = \/n — 2. 



n 


FIG. 5: The relations between n and H. (a), (b), (c), (d) 
and (e) were given by fixing up = 0.760,1.15,1.50, 1.80, 2.01 
respectively. Values the normal Friedmann equation holds are 
also indicated at n = 0. As the number of extra-dimensions 
increase, the expansion rate is suppressed. The end points of 
each lines are (n, H) = (n, 0.01). 


IV. CONCLUSION 

We have solved the Einstein equations and the scalar 
field equation for a global defect in higher dimensional 
spacetimes (> 6 ). The defect has a (3+1) dimensional 
core in n > 2 extra dimensions. We have extended 
the analysis by Cho and Vilenkin js] to other extra di¬ 
mensions and found numerical solutions without “cigar 
ansatz”. We have examined in detail the relation among 
the expansion rate H and the symmetry-breaking scale 
p and the number of extra dimensions n for these solu¬ 
tions. We find that even if the extra dimensions do not 
have a cigar geometry, the expansion rate H grows as 
p increases, which is opposite to what is needed for the 
recently proposed mechanism ; 6 | for solving the cosmo¬ 
logical constant problem. Finally we want to notice that 
our nondiverging solutions require fine-tuning of param¬ 
eters (n, up, (k/A 1 / 2 )#). So the problem of fine-tuning 
remain exist even if H decreases as 77 increases. 


fixed at some values. At n = 0 in this figure, the values 
read from the usual four-dimensional Friedmann equa¬ 
tion, H 2 = K 2 p/ 3, p = A? 7 4 /4 are also indicated. This 
figure shows that the expansion rate determined by the 
Friedmann equation is suppressed as the number of extra¬ 
dimensions increases and H vanished at a specific dimen¬ 
sion. This effect may be considered as the dilution of the 
cosmological constant. However, very small but non-zero 
expansion rate of H ~ 10 -33 eV cannot be reproduced 
without a fine-tuning for up. 
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APPENDIX A: NUMERICS 

In this Appendix, we give the details for finding the 
numerical solutions. 

1. (k??) 2 < n — 2 Case 

In this case, our strategy for the calculation is divided 
into 2 steps. 

1. Firstly numerically integrating Eq. ©, © and © 
as an initial value problem from the origin. 

2. Then solving them as a two-point value problem 
by the relaxation method with the initial solutions 
obtained by the previous step. 

In the Step © we solve the differential equations by 
the 4th-order Runge-Kutta method. The calculation 
starts from the origin with the initial conditions, B( 0) = 
C(0) = 1, B'{ 0) = C"(0) = 0,<£(0) = 0 and differ¬ 
ent 0'(O). The calculation is so sensitive to the initial 
condition, 0'(O), that bad choice would lead to the di¬ 
vergence of B,C, (j) at finite distance. So we have fine 
tuned the sixth initial condition 4>'(0) so that the point 
of divergence goes as far away as possible which we call 
rii m . Then we can obtain an approximate solution in the 
range [ni m , r m ax] without numerical divergence by fixing 
0(r) = 0(?’ii m ) at r > ru m artificially. In the Step© we 
solve the differential equations by the relaxation method 
with two point boundary conditions at r = 0 and at r max . 
The sixth condition is now replaced with f(r max ) = 0. 
After the iteration converges (typically relative error be¬ 
low 5 x 10 -10 ), we finally obtain a solution. 

2. (k??) 2 > n — 2 Case 

For (ft??) 2 > n— 2, arbitrary H including H = 0 leads to 
a divergence in B or C at finite distance from the origin 
and the singularity is formed. We call this point ?' S mg. 
At the specific value of H , the distance to the singularity 
becomes as far as possible and the divergence vanishes. 
We call this point rf. 

We solve the Einstein equations and the equation of 
motion of <f> numerically in the range of [0, r s i ng ]. But the 
value of r s in g remains unknown to be decided by solving 
accurately. So, our strategy of the calculation is divided 
into 4 steps. 

1. Firstly numerically integrating Eq. ©, © and © 
as an initial value problem from the origin. 

2. Then solving them as a two-point value problem 
by the relaxation method with the initial solutions 
obtained by the previous step. 

3. Extending the solution of (f> to sufficiently large 
value of r. 



FIG. 6: This graphs shows solutions which have singular¬ 
ity. The solid line and the broken line is for the parameter 
(n, «??, (k/X 1/2 )H) = (3,1.09,0.00375), (3,1.09,0.00400), re¬ 
spectively. These solutions have singularity at finite r. 


4. Numerically integrating the equations with the 
fixed 4> obtained in the Step© to find the true po¬ 
sition of the singularity and the solution. 

Step ©is executed in the same way as Step 1 of (ft??) 2 < 
n-2 case. But the range of r is [0, r s i ng ], where the bar 
means that this position of the singularity is different 
from the true one. Step ©is also executed in the same 
way as ( nrf) 2 < n — 2 case. But the boundary condition 
is imposed at r = 0 and at r = f s j ng — Ar. The right 
point is shifted to left by Ar so that the singularity is not 
included in this range [0,r s i ng — Ar]. By executing the 
Step© solutions of B , C , <p can be obtained and (f>'(r s i ng — 
Ar) is vanishing. So we can extend 0(r) to sufficiently 
large value of r, this procedure is the Step © In the 
Step© we solve the differential equations by treating the 
B{r ) and C(r) as unknown functions and <f>(r) as fixed 
background obtained in Step© The algorithm used is the 
same as Step © From the Step © to © we finally obtain 
the numerical solutions and the point of the singularity. 

Next, we find the ’’eigenvalue” under the condition 
that the position of the singularity becomes as far as 
possible. We can find the unique H as the eigenvalue 
with (n, ft??) fixed like as follows. If H is smaller than 
the eigenvalue, the energy of the defect core bends and 
closes the bulk space such that B(r) —» 0 as r —> r S m g , 
which is shown as solid lines in Fig. © Then C(r) diverge 
at the same point and forming the singularity. As H is 
increased, the energy density of the defect is consumed 
to inflate the brane and the bulk’s curvature is relaxed. 
If H is beyond the specific value, then that B{r) —> oo 
at finite distance r s i ng , which is shown as broken lines in 
Fig.© At a very specific value of H which exists between 
these two values, neither B(r) nor C(r) diverges and the 
position of the singularity has a local peak here. This 
H is to be called the ’’eigenvalue”. Typical examples are 
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(B2) 


M C (f ma x)ri 


2 

max 


(n — l)(n + 5 )(k?j) 2 
2(n + 2 )[(k? 7 ) 2 — (n— 2)]' 


(B3) 


where 


fc = 


I n + 2 (k??) 2 — (n — 2) 

2(n + 5) 2 nr] 


(B4) 


FIG. 7: The relation between r/ and H for cigar type solutions 
with n = 3. 


given in Fig. 121 and Fig. 0 


APPENDIX B: CIGAR TYPE SOLUTIONS 


In the previous works H, fg], the cigar type solutions 
in n = 3 case are studied. These solutions have asymp¬ 
totic forms y/Xi)Cr —> constant. The analytic solution 
for arbitrary n > 3 is 


ft 4 * ( r max) — 


2 (r 


4) — (n — 1)(kt/) 1 2 3 4 
(n + 5) 


(Bl) 


and r max is a sufficiently large value. 

It is concluded in jg] that eigenvalues in 
(kt), (k/A^ 2 )H) space can be line fitted. But we 
have reexamined this calculation and found more 
complex structure shown in Fig. [3 

Our method of the numerical integration is as fol¬ 
lows. We have set the seven boundary conditions B( 0) = 
0(0) = 1, B'( 0) = C'(0) = 0, m = 0, 0'(r max ) = 
0, (UnaxO(r max ))' = 0, and considered the equations 
©,©, © and H'(r) = 0, in which H is treated as 
the dependent variable of r. The strategy to solve these 
equations is the same as described in Appendix^ The 
iteration is not converged for some values of nrj corre¬ 
sponding to the blanks in Fig. 0 
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